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Introduction. 

It is well known that, if a rigid body whose principal moments of inertia are ^, ^, 
be set rotating about its axis of symmetry, and then be subjected to a slight dis- 
turbance, it will execute oscillations about its mean position, in consequence of which 
the axis of rotation will undergo periodic displacements relatively to the body in a 
period which bears to the period of rotation the ratio ^ : €1^ — ^. The object of the 
present investigation is to determine to what extent this period will be modified if the 
body, instead of being perfectly rigid, is capable of elastic deformations. 

The problem has important bearings in connection with the theory of the Earth^s 
rotation. The remarkable researches of Dr. S. 0. Chandler, published in a series of 
papers in the ' Astronomical Journal,''"^ have placed it almost beyond a doubt that the 
axis of rotation of the Earth is subject to a series of displacements, the most impor- 
tant of which consists of a periodic motion in all respects similar in character to the 
oscillation mentioned above, but having a period considerably in excess of that which 
theory would require if the Earth could be regarded as perfectly rigid. It is natural 
to suppose that this motion has its origin in the same cause, but that the theory by 
which the period has previously been assigned is in some respects defective. The 
most plausible attempt which has yet been made to correct this theory is that given 
by NEWCOMB,t who shows, by an elegant geometrical method, that the elasticity of 
the solid portions of the Earth and the mobility of the ocean will each have the effect 
of prolonging the period. He then proceeds to obtain a numerical estimate of this 
extension, basing his calculations on certain results given by Thomson and TaitJ 
with reference to the deformation of an elastic sphere. In order to make these 
results applicable, several assumptions have to be made which do not appear to mo 

* Vol. 11, et seq. For a summary of Chandler's resuUs, vide ' Science/ May 3, 1895. 
t * Monthly Hotices of the Royal Astronomical Society/ March, 1892. 
t ^ Natural PMlosophy/ Part IL, §837. 
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to be well founded ; and it is with a view to examining these assumptions that I have 
attempted to exhibit the sohition of the problem in an analytical form. 

The analysis in the present paper is confined to the case of a homogeneous spheroid 
of revolution, composed of isotropic, incompressible, gravitating material, while no 
account is taken of the surface waters. For further simplification we have also sup- 
posed that the figure conforms to that required for hydrostatic equilibrium, so that 
when the body is undisturbed we may suppose it free from strain in its interior. We 
have reason to suppose that this condition is approximately realized in the case of the 
Earth. 

In §§ 1-2, I have obtained the rigorous dynamical equations for the oscillations of 
such a system. The form of these equations is, however, such as to render an 
approximate solution necessary. Hence it has been assumed in the subsequent 
sections that the ellipticity of the spheroid and, consequently, the angular velocity 
of rotation, are small quantities. This assumption leads to a considerable reduction 
in the diflferential equations which express the elastic displacements, and, in fact^ 
reduces them to the familiar equations for the equilibrium of a strained elastic body 
The boundary conditions are also similar in form to the surface equations which are 
obtained in treating of the problem of the deformations of an elastic sphere, and as 
the solution of this problem is well known, we are in a position to obtain a solution of 
our differential equations applicable to the problem in hand. §§ 3-4 are devoted to 
the transformation of the equations into a form convenient for solution, while the 
actual solution is given in § 5. 

The raethod of approximation followed up to this point fails to lead to a determi- 
nation of the period. When, however, the body is supposed perfectly rigid, we are 
able to determine the period accurately by means of the equations of angular 
momentum for the whole system. This suggests the use of an analogous process, 
which is employed in § 6, to determine the period when elastic deformations are 
taken into account. 

The principal results of this paper will be found in § 8, where they are compared 
with the hypotheses made by Newcomb. It is found that the general character of 
the motion agrees with that assumed by Newcomb, but that his quantitative law as 
to the displacement of the pole, due to elastic distortion, is slightly in error. The 
bearing of these results on the theory of the Earth^s constitution is discussed in the 
final section (§9). 



§ 1. Differential Equations of Motion of Isotropic. Incompressible^ Mastic Solicit 

referred to axes rotating uniformly. 

Take as axis of z the axis of rotation, and let the system be rotating with angular 
velocity w about this axis. Let x^ + t^, y^ + t', % + '^^ be the coordinates at the time 
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t of tbe particle, which, when the body is -unstrained, is at the point % y^^ z^. The 
velocity components of this particle will be 

ft • 

U = — (a?o + t^) — cy {y^ ^v) = u — o^ {tj^ + v), 

Y = — (i/q 4" ^) + ^ (^0 Hh ^^) = f; + G) (a^o + te), 
W = — (% + ti)) = w. 

The components of acceleration will be 

* »»•• •»• 

If P, Q, R, S, T, U denote the six components of stress at the point x^ y^ z (where 
a^ = oJq + '^*? &<^«)? Xj Y, Z the components of bodily force at this point, and p the 
density of the material, the equations of motion may be written down in the 
same manner as if the axes were fixed, provided that we replace the accelerations 

dhi/df^ (Pv/df^ -dhDJdt^ by the valnes we have found above. Thus we have"^ 

^ «|, ^^ ^^ p^ ='p (14. — 2mv — mhi .— m^Xn)^ 

iJQO C/y G^ 

_- -|„ ™™- ^. — «^ pY = p (t» 4" 2mu — ct>% ' — ^ ft^Vo)^ 

3T 3S 0E , r/ 

ox dy az ' ' 

We have here employed two different sets of independent variables. The qnaa titles 
involved on the left have been supposed to be expressed as functions of the variables 
a;, f , 0, 1, that is to say, the coordinates of a definite point on which w^e fix our attention 
and the time ; these variables are analogous to the Etderian system in Hydrodynamics. 
On the other hand, the quantities % v, %v on the right have been regarded as functions 
of a?oi 3/0J %> ^j which correspond to the Lagrangian system of independent variables 
in Hydrodynamics. It is desirable for ns to retain only one set of variables ; we propose 
to select the former set and proceed to examine the modified form of the equations of 
small motion. 

If the symbol djdt be used as above to denote partial differentiation with respect 

^ LoYSj * Elasticity,* yoL 1, p. 60. 

MDCCOXCVI.— A, .2 T 
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to the time, on the supposition that x^^ y^, Zq^ remain constantj and 3/9/5 be used when 
Xy y, z, are the other independent variables^ we have 

^ __ ^ du 3 dv d ddo 3 
d^t dt dt dx dt 3;// dt 3-<;* 



Hence, if we neglect squares and products of the small quantities u, v, w, we 
obtain 



u 



dii/cU = du/dt, u -^^ cPujcW =: d^u/dt^, &c. 



Thus, the oaly modification necessary on the right will be to replace x^ -}- uhj x, 
and y()-\- V by y, and the equations of motion become 



9P , 3U 3T ^ 

+ gT + ;^^ + P^ 



del]' 

3U 



dz 
3S 



9^ + 3^ + S" + ^^ 






dx + ay + a; + ^^ 



pw 



^ 



(1). 



We have retained the fluxional notation, it being understood that the dots now 
denote differentiation with regard to the time on the supposition that x^ y, % remain 
constant. 

Since the stress-strain relations do not involve differentiation with regard to the 
time, they may be written down in the same manner as when the axes are fixed. To 
the degree of approximation to which we are going, we may replace dujdxQ, &c., by 
du/dx, &c., and therefore when the material is isotropic and incompressible, the 
components of stress are given by 



^ da 



Q :^ »~JD + 2!t 



tAJf 



dv 
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^ + 2n 'r.~~ 

^ ox 



R :::--— p -f 



3^0 

3^0 



^ /dw , dv , 



T ==n 



\dijQ ' dzr 

die dw 
+ 



It 



— p + 21t 

— • jp ~|- Z 

dv 



dv 
dy 

dio 
dz 



010 

dy 



■BX&RUHl 



dz, 







dx 



0; 



V 



u 






ov 



— p 



a 



u 



\dx^ 83/, 



dz 

dio 
dx 

dv dii 
dx dy 



Idib 

tt j ,;. ~ "4- 

oz 



J 



(2), 



where p denotes the hydrostatic pressure at x^ y, z, and It the rigidity. 
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QOQ 



If fortlier, the bodily forces be deriYable from a potential function V, so that 

X = BY/dx, Y = BY/dy, , Z = BY/Bz ,..•.. (3). 
on substituting the values (2). (3), in (1)^ we obtain 



n 



B 






p J 

z 
z. 

p 



u 






= t; 4" 2(x>u^ 



w. 



Putting tt/p = % V + i«^ {^^ + ^^) "^ pIp ~ #? these equations take the form 






-^ 



« . * 



• » • ♦ l4r/» 



J 



We must in addition express the feet that the material is incompressible; this 
is done bj the equation 

BujBx + S'y/^l/ + 3^/9^ = , . . . . . . . (5), 

The equations (4), (5), ■which are the rigorous equations for the vibrations of a 
rotating^ incompressible^ elastic solid^ are theoretically sufficient to determine u, v^ tv,\P% 
subject to certain boundary conditions. Eliminating v^ u in turn from the first tw^o 
of equations (4), we obtain 



nV 



% 



nV 



s 



Bi^ 
Bt^ 



s 



+ im 



+ 4ct)^ 



8^^ 
Bt^ 



lb 



nV^ 






V 



nV 



% 



Bf J Bw 

B^ \ 3i|r 
Bfi J By 



• "-f- MO} rs. 






1, 

8^ 



Byj' 

8^ 



Applying the operators BjBx^ ^l^y^ ^nd making use of (5)^ 






82 N3 



ot 



0, 9^ 

*"T"" ^ ffj /%' r\ 



Bti) 
80 



n 



3 



03 \ /gs^ 8^-^' 



Bf 



\8a]^^ ^ 



8^'^ / * 



Hence^ on eliminating w by means of the third of (4)^ we obtain the following 
equation for ^ :- 

\j • • « • • • iij/» 



83 \B 94 



3^5 



8^^.8# 



t|i 



9 rp 9 
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Ill the case where n is zero, it will be noticed that the last equation reduces to 
Poikoarje's differential equation for the oscillations of a rotating mass of liquid.* 
It may easily be shown by retaining any one of the four quantities % v, iv, \jj^ and 
eliminating the other three, that each one of these quantities is a solution of the 
equation (6). 

If the motion relatively to the moving axes consists of a simple harmonic vibration 
in period 27r/X, we may suppose t^, v^ tu^ \}j each proportional to e^^^, and the equations 
(4), (5), (6) will then become 



•^ 



(n V^ + X^) u 4" 2miXv := — dtf^/dx, 
in V^ -}- X^) t; •— 2mXu = — 3^/3^? 
^^^ V^ + X^) tv := — dijj/dz^ 

dii/dx + dv/dy-l-dw/dz =:= 0^ 



in 



(8). 



I 2. Boundary Conditions. 

The cmiditions to be satisfied at the boundary are that the components of surface- 
traction should vanish at all points of the displaced surface. We proceed to replace 
these conditions by certain analytical conditions at the m.ean surface. 



.^■^ 






/ 



./ 



y 



y 



/ 



Take a point P on the mean surface and let the normal at P meet the surface of 
the distorted body in P^ 

Let cos a, cos ^^ cos y be the direction-cosines of the normal PP^ and cos a + \^ 
cos ' ^ + Wj, cos y 4" % the direction-cosines of the normal at P' to the displaced 
surface ; also let PP' = C- Then li, Wj^ %, C will be small quantities of the order of 
the displacements u^ v^ %o. 

Let P, Q, E, S, T, U denote the components of stress at P and the same letters 



* * Acta Mathematical vol. T^ p. 356. 
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accented the components of stress at P', and let cM denote an element of the normal 
to the mean surface. Then to our order of approxhnation we have 

P' =z p + ^ d?ldn\ Q' = Q +• c ^Q/^^/^; &c. 

Thus the component of surface- traction at P^ parallel to the axis oi x is 

P^ (cos a + ^i) + U^ (cos P + 'W^) + T^ (cos y + ^^i)j 

/ ?iP\ f r)TT\ / f(V\ 

=^\^ + ^ ^'^') (^^^ ^ + ^i) + (^U + I g-j (cos y8 + m^) + {T^l g^j (cos y + n^), 

/3P 3U 3T ^ 

= P cos a + IT cos y8 + T cos y + H ^/ cos a + ^, cos ^ + ^ cos y 

if we neglect small quantities of the second order. 

Now in the small terms we may replace P, Q, R^ &C.5 by their values in the 
steady motion. Since we have supposed the body when undisturbed to be free 
from tangential stress in its interior, we have in this case 

PzzzQzzzRzzz ~~_^9, S = Tz=:U=;0 

throughout, while at the surface also j9 = 0. 

Therefore IjP + mjU + ^iT, dXJ/dn^ dT/dn may be put equal to zero, and 
;-- _ _ =: _. p — - ^ Y -|. 1.^2 (^x^ ^ y^'j'j :=z ^ pg say, where g denotes the value 

of gravity (inclusive of centrifugal force) at the surface. Thus the ir-component 
of surface-traction at P' is 

P cos a + U cos i8 + T cos y + pgl cos a. 

Introducing the values of P, U, T in terms of the displacements and equating this 
expression to zero, we obtain 

— p cos a + tl -rT" cos a + -;r- cos + -;r~- COS y] + it hr" COS a + -^T" COS p + TT- COS y 

^ \dx ay ' dz ' ) ' \dx '3^ ^ ^ Zx ^ j 

= — pgl^ cos a ; 

and, in like manner, by considering the components of surface-traction in the 
directions of the axes of y and z^ 
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— p COS p + tt ( 'V- COS a + -X- COS p + ^ COS y 1 + It ( -^ cos a + -^ cos p + "^r^ cos y 

= ■- pgC cos /3, 

— p COS 7 + It v^ cos a + 7^— cos p + "^ cos y 1 + It [ ^ cos a ^- -^ cos p + aT ^^^ 7 

=:== — pg'^ COS ye 

These equations express the fact that the surface-tractions at the mean surface are 
equivalent to a normal stress equal to the vraght of the harmonic inequality, 
and they might have been written down at once from this consideration. It has, 
however, been thought preferable to verify them at length ; it will be seen that the 
shorter procedure is only justifiable in the case w^here the material is initially in 
a state of hydrostatic equilibrium. If in the zero configuration this condition is not 
satisfied the form of the boundary equations will be much more complicated. 

Let us now replace p) by the function \jj of the previous section. By the definition 
of xjj we have 

pIp zz: V + i <^^ (^^ 4- y^) "- ^P 

=== non-periodic terms -{- v^ — xjj, 

where v^ denotes the potential due to the harmonic inequalities. 

The non-periodic terms vanish at the surface in virtue of the conditions for steady 
motion, and therefore the boundary equations may be written 



xjj cos a + 7^^ 



TT- COS a-+- TT- COS p + ^ - COS y + ;:r" cos a + V cos p + K- COS y 

ox cy oz ' ox ox ox ' 



> (9), 



~: (v^ — gt) COS a 
xjj COS p -j- n i^ COS a + -r.^ cos p -r -^^ cos y + ^ cos a -^ -^ cos p + -^-" cos y 

— (^' — gl) cos ^ 

, , \dw . d'W ^ , dio du , dv ^ , dio 

xp cos y + '^ ]^ cos 06 + Y^ cos p + v" cos y + ^ cos a + ^- cos p + ^ cos y 

:= ('y' -— gQ cos y 

where 

^ = '^^ cos a + 'y cos ^ + ^(; cosy (lO)* 

The rigorous method of procedure would be first to solve the equation (8) for \p. 
Having found \p we could introduce its value into the right-hand members of (7) and 
proceed to find solutions of these equations consistent with the boundary conditions 
(9). Unfortunately the form of equation (8) is such as to make it appear hopeless 
to carry out this process, and we must have recourse to some method of approxima- 
tion before we can advance further. 
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§ 8. Change of Variables, 

An elastic body, sucli as that with, whicii we are dealing, will, of course, be capable 
of an infinite number of independent normal types of vibration. _ The equations of 
motion we liave found in § 1 are applicable to any one of these types, while the 
substitution of the solutions in the boundary equations should lead to an equation 
for the determination of the frequencies. In general, the values of u, v, w will 
become very small when n is large in such a manner that nu, nv, mv approach finite 
limits, while the admissible values of X become large of the order n^. That type of 
oscillation wath which we are concerned is, however, unique in character in that it 
continues to exist even when the rigidity is perfect. If^ then, in the expressions for 
u, V, w in terms of n we suppose n to be made infinite, u, v, iv should approach finite 
limits which we will denote by ii^, v^, iv^. 

The quantities ?^q, v^, Wq denote the displacements of a body w^hich is supposed 
perfectly rigid. Now the most general small displacement of such a body consists of 
a translation whose components we may denote by ^, rj, ^, and a rotation whose 
components we denote by d^, 6^^ 6^. Thus, the most general values of Uq, Vq, w^^ 
geometrically possible are 

When, however, we are deahng with the rotation of a rigid body not subject to 
external disturbing force the quantities ^, rj, L ^3 will not appear, and we may take 

Uq = z6^, Vq = — zdi, Wq = — x6c^ + y^i • (ll)- 

Let us now suppose that 

qi zz: iCq + Uiy V := Vq -{- V^^ IV = 10^^ + W^, 

where Uq, % iVq have the values (llj. This is equivalent to supposing that the body 
is first displaced by rotation as a whole through small angles d^, 0.2 about Ox, Oy and 
is then subjected to elastic distortion, the displacements due to the distortions 
being 11^, v^, w-^. 

Further, let C = ^0 + ^u '^^ ==^ '^'o + '^'u w^here 



^Q = qi^ COS a + t?Q COS /3 + Wq cos y 
^^ = it^ COS a + ^1 cos jS + '^^1 cos y 



■^ 



r" . . . . ' . . ( J. -U ), 



rJ 



and where Vq, v\ denote the parts of v due to the harmonic inequalities ^q, ^^ 
respectively. 

If Vq denote the non-periodic part of V, Vq + Vq will be the potential at the point 
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Xy y, z due to the attraction of the body when it is rotated without distortion through 
small angles ^i, ^g* If then x^, y^^ z^ denote the coordinates of this same point referred 
to axes obtained by rotating the old axes with the body, on putting in evidence the 
arguments of the function Vq, we have 



Vo {x, y, z) + v^ = Yo {x^, 3/1, z{} 



» f » e • 



• (13)- 



But the direction cosines of the two sets of axes are given by the scheme 





X 


y 



z 


tJb ] 

2/1 


1 




1 


1 


h 


0, 



whence x^^ :=z x — z0o, 2/i = 2/ + ^^v ^1 = ^ + oc0^ — yO^. 

Introducing these values in (12) and expanding by Taylok's theorem^ we find 



"^0 (^^ y^ ^) + '^'o = 'Vq{x — zO^, y + zO^, z + xd^'- yd^) 

avn . . av 

dy ' ^ '^ " ^' az 



= Vo {x, y, z) - ze, -^ + zO, ^ + {xd, - ye,) ^ 



But from the definition of g we have, at the surface. 



g cos a 



d\ 



r 







0) 



IX 



%', g cos /3 :^ 



av^ 







^y 



oj-^y^ g cos y 



and, therefore. 



ay, 



v'(^ = Z02 (g cos a + o)^x) "^ 2^1 (g COS ^ + o/y) — (xO^ — y6^). g cos y 

= ^ (t^Q cos OC -^ Vq cos ^ + '^o <^^^ y) + ^^ (^3^'^ "~" ^iy^)i 



or 



whence, finally, 






f e 9 * 



• (14) 
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If now we change our variables from u^ v, iv to iii, Vi, lOi, the equations of motion 
become 



nVhti + \hii + Sco^Xt^j = — d^/dx •— X^^^ + 2mkzdj^ 



~\ 



nV^i 4" ^^'^1 "^ 2mXi(n 



dxjj/dy + Xh0i + 2m\z0.2 



nVht\ + ^^^^^1 = "^ d-^jdz + X^ {006^ "^ V^i) 



1^ • • (15), 



->' 



while in virtue of (14) the boundary equations may be written 



du. 



ti/ COS a 4- n i ir-~ cos a + -^ - cos p + -^ cos y + :;;— cos a + ^ cos 8 4- -^ cos y 

^ ' d:^ d?/ & 'da? d^ ' d^^ ^ 



— {v\ -— ^g^i) cos a = 0)^ {0^^^ — ^i^^) ^^^ ^? 

ti cos p 4- ^ 1 ^r" cos a 4- ;. - cos p 4- x~ cos y 4- "^ cos a 4* ^ cos B 4 -^ cos y 

^ ' [da3 ' oj/ ' ' d^ % f^|/ ay ^ 

— ('?/| — ^'Ci) COS fi =: 0>^ {0%^^ — ^i2/^) cos ^, 



li COS y A- n ^ -^ COS a 4- "^ COS p 4- '^ COS y + tt^ cos a 4* -^ cos « 4 —^ cos v 

^ ' ox otf oz ^ oz az oz ^ 



- (^'1 - 9C1) cos 7 = 0)^ {0^xz --- 6^i2/^) cos y. 



>- 



(16). 



^ 



I 4. Rechiction of the Equations tvhen the Body is of the form of a Spheroid of 

Small Ellipticity. 



The only assumption we have made as yet as to the form of the free surface is that 
it is a possible figure of equilibrium for a rotating mass of liquid, so that the body 
may be free from strain when rotating uniformly. The simplest form which can 
occur and that which presents the greatest interest is the case of a spheroid of 
revolution of small ellipticity c. We propose for the future to confine ourselves 
to this case. If we neglect the square of € the angular velocity of rotation is related 
to e by the equation 

» • • • » • • • * iiffj 



€ = 15ct)7l67rp . 



where the density p is expressed in gravitational units. If then the units of length 
and time be so chosen that p is finite, m will be a small quantity of the order el 
Take as the equation to the free surface r = a { I 4" ^TpJ where 



MDCCOXCVI.— A, 



T.^ = {x^ + y^ --^ 2e^)/3a^ 

2 u 
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Then the direction cosines of the normal to the surface t= a {1 + ^Q/*}? where Q., 
is a solid harmonic of order n^ are 



1 



J . I fl/d. 






> 



J OwO. J 



and thus we have 



cos a 



cos ^ 



cosy 



m 
r 


+ ae^ 


V 

r 


-- ae < 


z 
r 


+ ae - 



Zi'jj 






%r 



ax J 



a -^3 






3t„ 

3T( 



whence 



c 







^.^3 cos a — 2;^;^ cos ^ — (O^x — ^j3/) cos y 



ae 



ae 



'h 


\ 


3T, 


- z 


9.r y 


)+ 


«-. 


^. 


re" 




«,. 


2^^ 



a-' 



dy 



z -^ — -^-"^ 2/ 



9T 
9s; 



or 



U=-lw.-eM 



a « • \ i- k) 1 1 



Now observation indicates that in the case of the Earth the oscillation in question 
differs but slightly in type from the motion of a rigid body, whence we conclude that 
u^, Vi, tc^i must be small compared with u, v, m We propose therefore to make the 
assumptions, leaving the verification thereof to our subsequent work, that u^, Vj^, lOj, xfi 
contain the small fector m% while X is a small quantity of the order o)^. The latter 
assumption is justifiable when the body is perfectly rigid, since in this case we have 



€0 






If then, we neoiect small quantities of the order m^ or e^ the equations of motion 
(15) reduce to 



nVhi, 



9i|r 



d-JG 






JL — i J ■-- 

9^ 9^ 



9l|r 

= 



nV^w. 



3^ 



1Z 



L . . (19). 






Since all the terms in the boundary equations involve the small factor m^, we may 
replace cos a, cos A cos 7, by .^/r, yjr, zjr respectively, with errors only of the same 
order 0)^ and thus the approximate form of the boundary equations is 
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oo X 



,. . ,. (. ■ 



9_ 

3?' 



d 



= X. ai^ (d^jjcz — ^jj/z), 



«--.9Ci)a? 



y ■ w^ (^2a:;2 - 0^y7), 



z . ft)^ (^gccz — - 6iyz). 



^ 



/ 3 \ 3 



> 



/ 1 \ r) 



J 



(20).^ 



To the same order we may suppose these equations to hold good at the surface r = a, 
mstead of at the surface r = a {1 -f ^T2}- 



§ 5. Determinatiofi of the Elastic Distortions. 

The modified forms (19) of the equations of motion are dimply the equations to 
which we are led in determining the displacements in a strained elastic solid and the 
solutions of them, when either the displacements or the surface-tractions at the 
surface of a sphere are given, are well known. t We may readily adapt these 
solutions so as to satisfy the boundary conditions (20). 

Denote for brevity the function co^ [d^xz — Oiyz) by Sg. 

From (19) we obtain at once V^i/; = 0. This equation replaces the more com- 
plicated form (8) ; a particular solution of it is i/i = AS^. 

Introducing this value oi xjj into the left-hand members of (19), we find 



V%1 



A as. 

71 d-JG 



V%. 



A dSo 

ji 

n dy 



, Vhi\ = 



A 3Sc 



n 



3; 



particular solutions of which are 



u. 



10 



A o 3So 
n ' dx 



i'i 



10 



n oy 



W, 



10 



/v*."^ it 

n dz 



• (21)' 



These do not satisfy the last of equations (19), but they make 



dx dy dz 



2 
5 



A 



n 



a 



0* 



* For the physical significance of these equations, vide § 8 infi^a. 
t Vide LoYE, ' Elasticity,' vol. 1, cliap. 10. 
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We must therefore add to the particular integrals just found complementary- 
functions which satisfy the equations 



VV.1 = 0, V\ = 0, Vhi\ = 

A 
3^ ' di/ ^ dz n 



"^ 



+ 1 5 s^ 



> 



a • » « • • \ Zj Zj j t 



Now the functions r'^ tt- ( -^ ), ^^ ^ (^ )? ^^ ]v" (7;^ ) are spherical harmonic func- 



tions which remain finite at the origin, and hencOj if we take 



%u = 



'd. =^- 



B i^ + Cr^ -^ 
d=x? dx 

di/ dy 



9 /S.. 



T>!) 



tfi 



oz d, 






V 



1^^ 



• (23) 



^ 



the first three of equations (22) will be satisfied identically, and the fourth will also 
be satisfied if 



3.7.0 



5 



A 



n 



or 






10 5 






n 



Therefore, adding together the pai'ticular integrals (21) and the complementary 
unctions (23), we obtain as a set of solutions of the equations of motion which remain 
finite at the origin 



I" 



t|; = AS3, 



'-S 



U- 



V, 



^6'i = 



1 r: n^^ ^3 



10 



71 



. a. 



G 



+ B 



-^■2 






as. 



10 



10" 



71 'dy ^ "^ 3^ 



A .3Sc 



'?' ''?r~ -p XJ 

71 OZ 



y{ 



_^i_ A ^..7 z^ 1^1-- 

10 5.. 2^ I ^.5 



10 5 



A 






Sc 



«-»•) 



A „ 3 /S. 



3.,^ 



X VJ «5 



?^ 



??- 



3,i! \ T^ 



> 



■ ' (24). 



ii 

J 



We must now prove that these solutions are suGficiently general to satisfy the 
boundary conditions (20) at the surfece of the sphere r ^ a. 
From (24) we obtain 



u^x + "^iV + ^h^ 



X A ,.2s j^ 2m, + fs ^ r^S, 



A 



n 



n 



1 - r^Sa + 2BS2 . (25). 
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Since in the small terms which already contain S^ as a factor we may treat the 
spheroid as a sphere of radius a, we have 






2B - i ^ • 

n 



s. 



and therefore 



V- 



I Trpa^i = I TT/) 



2B 



1 

7 



ka^' 



n 



Soj 



whence 



^i' — 9^ = (s ^/>^ — 3 ^/>^) ^1 



-h ^P 



2B-i 



Aa^-i 



n 



s. 



(26), 



(27). 



Consider now the different terms of the left-hand members of the boundary 
equations (20). By means of the formula 



.^2 7\^ 



JUy^O "—— ^ r\ ■"" 



we obtain at the boundary 



X\JJ ZZ2 A 



5 Ox 



cc 



5 8.:^ \ r^ 



asg _a7 3 /vj 



5 dx 



a 



o dx \7 



^0 



Also we have 



^'a;. - 1)^1 



, A o ^^3 



n 



da 



„4_ 
105 



J ' 



A . 8 /S. 



7" 



7^ 



8^ \ r 



5~ / ' 



and from (25) 



9^: (^1^ + ^1?/ + ^^1^) 



^ ^ ' n ox ox 



2 
7 



therefore, at the surface. 



A 



n 



'2 
r3 8S 



2 



5 8 



X 



5 8<2? V r^ 



J 7^ 0.?? 



I («,x + ^;32/ + ^.,.) = - i -^' ^^ + 2B ^^ + A 



A „ a /s 



a' 



■K. 



a 



X\ T 



dx ' 



Lastly, from (27), 



^ (^/ - ^'^i) = - 7-5 ^P 



2B-^ 



Aci/' 



n 



a' —' - a7 



S\ *n 
2 



L. 



8 



X 



dx \ r^ 



and therefore the first boundary equation becomes 



dx 



a'' 



A ^ - iAa^ - iAa^ + 27^B + ^Trpa^ 

8 /S. 



o 



7 



Aa3^ 



^ 



+ a^ 



?X' \ ^ 



-.5 



A-^ 



a" 



x^-Aa-^ + -3\Aa2 - -^V^a^ 



2B 



0.^ asg 

5 dx 



di _a_ / Sf,^ 

5 dx \ r^ 



1 A«2^v 
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The remaining boundary equations can be written down by replacing x hj y, z 
respectively. Hence, all three boundary equations will be satisfied provided A, B 
are subject to the relations 



- Aa^ + lOnB + -^^wpaN 2B 



JL 

7 



A.a' 



n 







Aa^ - /f Aa^ + i\rrpa^ ( 2B 



1 ^^'\ ^2 



n ! 



Solving these equations, we obtain 



A "^^^^^ ~" 



21 



19 + l'^ef^ 



3 ' 



Bn 



^a? 



19 + A^^fi^: 



O J 



n 



2B 



7 



Ati^' 



ba^ln 



n 



19 -f 4 



TTpa 



n 



Now, let e' = 5co^a^/38^^, so that e' will denote the elHpticity which would be 
induced in a sphere of radius a by centrifugal force when distortion is resisted by 
elasticity alone,* Then the above values may be written 



XJL 






19 {1 + e'ie}' 



Oi' 



8e75 



// „ J 



1 -I- e/e 



/ o"R 



Ao?\ 



\ 



n 



9J 



0) 



1 4- e'/e 



Finally from (24) 



u ■=: Uq -{- u^ 



i I V 42 3 ^ I 21 2 



.0 "S 






1 + e'/e 



^ "^ ^ 1 + e7e 



e' X (6^xz — ^1^;^) 

3 



a' 



V z=z Vq + Vi =: — ;^^i < 1 + _^ 



ie' 










1 + e'/e 



, 4 „_1 ^i^^ 

1. "T" € JG 



e^yz) 







'^ = u?o + i(^i = (•— 0:^2 4" 2/^i) ^ 1 









r e 



"N 



X "T~ €. J€ 



y (28), 



^ ^ ^ 1 + e'le 



a? 



^ 



J 



' * 0/ Thomson and Tait, "Natural Pbilosopliy/' Part II., § 837. 
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and from (IS), (26) 

C == Co + Ci = - 7 {0,ooz - 0,yz) + ~j -^':^^^ {0^ - 0,yz) 



26//, ^ V . 2 € 



2 



J ™ 

a [ 1 +- e /e 



(l7ot(/2' """~ ifj ) • * * • * ^ * ' \iUt/ J* 



S 6. Determination of the Period. 

The method we have followed hitherto has enabled us to express the displacements 
at any point of the body by means of two arbitrary constants 9i^ 6^^ ^^^ ^^® 
quantity X, whose value it is our chief object to determine, has entirely disappeared. 
We have, in fact, verified that the equations (7), (8), and the bouodary conditions (9) 
are approximately satisfied by the forms (28), while to the same order zero is the 
approximate value of X. We require now to have recourse to a method which will 
enable us to carry our approximations to the value of X further. 

By using the well-known equations of motion of a body of changing form Professor 
Woodv^ard'^ has shown that the determination of the period may be reduced to the 
evaluation of the disturbing angular velocities due to the flow of the material 
relatively to the principal axes of the body. As we have now expressed the 
displacements, and consequently the velocities, at any point in terms of 0^, ^3, 
wdiich may be taken as the displacements of the body as a whole, we ai^e in a 
position to calculate these disturbing angular velocities. We might then introduce 
their values in Professor Woodward's equations and pro ceed to the determination 
of the period by his method. We propose, however, to make use of equivalent 
equations which express that the rates of change of angular momentum for the 
system as a whole a-bout the axes Ox^ Oy are zero. It seems somewhat preferable, 
on account of the additional simplicity of the motion of the axes themselves, to 
refer to these axes rather than to the moving axes used by Woodward, which 
correspond with our axes Ox^^, O^/^. 

If hi, h^ denote the components of angular momentum about Ox, Oy, we have 



h 



{toy — ('y 4. xm) z] dm, 



where dm denotes an element of mass, and the integral is taken throughout the 
volume contained by the displaced surface. Replacing the integral by an integral 
taken throughout the mean volume and a surface integral, we have 

{(^loy — vz) — mxz] dm — p ^coxzdS, 

^ ' Astronomical Journal,' xv., No. 345. 
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or since 



Similarly 



X7^ dm \x% dm = and v, iv, t, are each proportional to e 

I (ivy — vz) dm — po) tpcz cZS^ 



hn 



X^ {ivy — vz) dm — pmiX tpcz dB. 



h^ = iX (uz — • tvx) dm -— po) | [,yz dS, 
J J j J J 



h = 



X 



3 



(ivy —vz) dm — » pwiX iyz dB. 



The equations of angular momentum are 

« 



]l^ -j- Jl^ff^ ^z 0. 



• • 



Replacing h^, h^, hi, h^ by the vahies we have just found we obtain the following 
rigorous equations 



2 



- X^ 



(wy — vz) dm — pmik tpoz cIS 



— iko) 



{i(z — ivx) dm + po)^ tjxz c?S = 



V a 



f* /• 



(t^^ — wx) dm — ■ pcoeX ^^2j cZS 



^ 



r r r re 

+ tXoo (tt??/ — 1)2;) (im — po}^ ^xz dB = 0, 

J J J J » 



» « 



(30). 



Now by using the approximate forms (28) and denoting by M the mass of the 
spheroidj we have 

r r f r r r 

(ivy — vz) dm ■■= { — - xy6<^ + (^/^ + z^) 6i} dm + {wiy -— Viz) dm 



J J 



2 



Ma^'Oi + terms of order e. 



{tiz — wx) dm = {(2^^ + *'^^) ^3 "^ ^2/^1} ^^^ + i^h^ "^ '^^'i''^) ^^^^ 



= f Ma^^o + terms of order c 



Also from (29). 



t,yz dB = 



2^/ e^ 



yh^ dB = ^iBTT-a^Oi i 



\1 + eVe/ ' 



^0:2; (iS = — -f-^7ra^d.2 



//- ^ 



1 + e'/e 



tlie errors being of the order e^. 
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Hence neglecting terms of order m^ only (X being regarded as of order co^) the 
equations (30) reduce to 






? 



+ iXco . ^Ma% + po^W. . '^\na^ - -''-77- = 0. 

1 -f- e /€ 



Tiiese equations will be consistent if 9^ = ido^ 



ft .^0 ^ 

-fr,7rpa''a}'' q — -—77 

|,Ma3«) ^"^f+^/e' ...... (^1). 



This gives the value of X with errors of the order co^. 



§ 7. Numerical Values. 

If we suppose the rigidity of our body to become perfect we should obtain e' = 0, 
and therefore X/co = e. 

This is the value we should arrive at if we started by neglecting the elastic distor- 
tions. We see now that it is too large and that consequently the effect of elastic 
deformation is to diminish the frequency or to prolong the period. 

The expressions for €3 e are 



iDTTp oSn 

Taking the sidereal day as 86164 mean solar seconds and using Boys's values'^ for 
the mean density of the Earth and the constant of gravitation, viz. : 5*5270 and 
6 '6 5 76 X 10~^ we find from the above formula that for a spheroid of the same mean 
density as the Earth, rotating in a sidereal day^ 



€ 






Again taking the Earth's mean radius as 6-371 X 10^ centimetres and 
tt= 8'19 X lO^^t which is the rigidity of steel, we find 

71— n/p = P482 X 10i\ 

^ ' Proc. Boy. Soc.,' 1894, p. 132. 

t Eyeeett, 'Units and Physical Constants/ pp. 61, 65. 
MDCCCXCVI. — A, 2 X 
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whence 



1 

522 



X/o) 



1 



"■ 1 + e'le ' 335 

We conclade that for a homogeneous spheroid of the same size and mean density 
as the Earth J the period would be extended from 232 days to 335 days in consequence 
of elastic distortions, if w^e suppose the rigidity to be that of steel. 

If we take into account the variations in the density of the Earth^s strata the 
problem, presented becomes much more complicated. We must replace e by the 
Precessional Consfcantj which will no longer be equal to the surface ellipticity. Its 
value may however be accurately determined by means of data furnished by the 
Theory of Precession ; this value is known to be 1/305, Hence the effective value of 
€ is diminished in the ratio 232 : 305. 

As regards the effect of heterogeneity on the value of e, we are, at present, only in 
a position to make speculations. Professor Newcomb points out that in calculating 
the mean density^ greater weight should be given to the density of the superficial 
layers on account of their greater effective inert ia^ and hence e^ should also be 
diminished. A reasonable hypothesis seems to be that it is diminished in the same 
ratio as e If w^e make this hypothesis, we find that if the effective rigidity of the 
Earth were as great as that of steely the period of the Eulerlan nutation would 
become 

— — days = 440 days. 

This period is slightly in excess of Chandler's observed period of 427 days. We 
therefore conclude that the effective rigidity of the Earth is slightly greater than 
that of steel. 

If we make the same hypothesis as above^ with regard to the effects of the 
variations of density, we may easily calculate what degree of rigidity would be 
consistent with Chandler's observ^ed period. We find for the period of a homo- 
geneous spheroid of the same degree of rigidity 

—--- — -days = 326 days. 

Putting X/m = fig, € = -^^-^ in (31) we obtain 

^ ~r^/^ -— ' "2 3 2 

e — 232^ ^ — 572? 
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and therefore 

^849x572^ 10" =8-98 X W\ 



§ 8. Physical Characteristics of the Motion, 
The height of the waves at the free surface is given by the formula (29), viz. : — • 

^ = - T r+7> (^^^^ " ^1^'^- 

Comparing this with equation (18) we see that ^ may be obtained from ^q by 

changing 6^ 9^ into 6^ 7, 9^ ;. Tiius, to our degree of approximation^ the free 

surface will remain a spheroid of revolution of elhpticity €. The position of the axis 

c 

of this spheroid may be found by rotating the axis Oz through small angles d-^ 



e + e" 



9^ — ^ — , about Ox, Oy respectively. 

Again^ from (28) we see that with a high degree of approximation the displace- 
ments at the point x, y^ z will be given by 

It = z9.2i V •=- — z9^, IV = -- x9c^ + y9^ (32), 

in other w^ords, the displacements due to elastic deformation will be negligible 
compared with the displacements due to the rotation of the body as a whole. We 
have seen, however, that the elastic distortions will produce an appreciable effect in 
displacing the axis of figure. 

The modified forms (32) indicate that 9^, 9^ are the displacements of an axis 
sensibly fixed in the earth, which we may call the mean axis of figure, while the 
motion, at any instant, will consist very approximately of a rotation, as a rigid body, 

• ft 

with angular velocity o) about an axis] whose direction-cosines are 9iJo}, 9Jo), 1. This 
axis we shall hereafter refer to as the instantaneous axis of rotation. It is, of course, 
only when we neglect the displacements due to elastic distortion that such an axis 
exists. 






We have found that when X = o) ( ;), 9^ = i^^. Taking 



where ^, r ai'e real, we have 

Zi X. Zi 
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Adding to these the solutions obtained by changing the sign of i^ wherever it 
occurs, we find as the real solution corresponding^ to the frequencies i o) -, 



^ -=. 2(j& cos X[t — t), 



0) 



26 -^ sin X{t — t), 

CO 



0. 



6^ ^=- 2<p sin \{t --- t), 

X 

w 



O) 



+ 2(f) - — cos X (t — t). 



y 




c 



M 



Let C be the point where the axis Oz cuts the surface of the spheroid j and let us 
take a pair of rectangular axes Gx, Cy, having C as origin and parallel to the original 
axes Ox, Oy, Let Pj P', R be the points in which the mean axis of figure, the axis 
of the deformed figure and the instantaneous axis meet the plane Cxy» Taking the 
radius of the spheroid as unity, the coordinates of P are 0o, — ^i, or 



2^ sin X (t — r), 



2(f} cos X (^ — r). 



The coordinates of P' are 0: 



^e -t- e^ ' 



1 e + e' ^ 



or 



2^ rfr- siii ^i^- ^)' 



2<p , cos X (t — t). 



and the coordinates of R are OJo)^ ^2/^? ^^ 



e3 



_ 2(h — ^ sin X (^ — r), + 2(^ -»_^— eos X {t — r). 



We thus see that the points P, P'', R all lie on the same straight line through C, 
and that this line revolves about C with unifortn angular velocity X relatively to the 
moving axes Cx, Oy. These axes are themselves rotating with angular velocity o) 
about C in the same direction, and, hence, the actual angular velocity of the line 
PP'"R about C is 0) { 1 + e^l{e + e')}. The distances of R, P, P' from C remain constant, 
while since OR is of the order e compared with GP or OP', which are themselves 
small quantities, we may suppose the point R to coincide with C. We conclude that 



MR. S. S. HOUGH 01^ THE ROTATION OF AN ELASTIC SPHEROID. 341 

the axis of rotation remains sensibly fixed in direction while the mean axis of figure 
and the axis of the distorted figm:'e describe cones of revolution about it in a period 
approximately equal to the period of rotation, these three axes always remaining in 
a plane. 

Let us next consider the motion of the instantaneous axis of rotation relatively to 
the Earth. Take a new set of axes O^j^j Oy^., Oz^ found by rotating the old set 
through small angles 0^^ 0^ about 0^, Oy ; these new axes will be sensibly fixed in 
the Earth. The motion at any instant consists of a rotation whose components about 

the old axes are 0^, 02, (o ; resolving these rotations about the new axes, we find for 
the components 

0^ — 0)0,2:, 02 4" <^^]J <^* 

Hence the direction-cosines of the axis of rotation relatively to the axes Ox^, Oyi, Oz^ 
are 

{0^-oj0.)/oj, {0^ + (o0^)/a>, 1, 



or 



'^ \ / 9 



2<^ ( 1 + rri'j '"' ^ (^ ~ "■)' +2c^(^l + -,^,-;)cosX(i-T), 1. 

Thus, relatively to the Earth, the axis of rotation will describe a cone of revolution 
about the mean axis of figure with angular velocity X, the direction of motion in 
this cone being the same as the direction of the Earth's rotation. This motion would 
manifest itself by a periodic change in the latitude of places on the Earth's surface, 
as found by astronomical observations. 

The circumstances we have here described are the well-known characteristics of the 
Eulerian nutation,"^ with the additional feature that the axis of figure is displaced by 
centrifugal force towards the axis of rotation. This latter fact has been assumed by 
NEWCOMBt and other writers;|: as the basis of their work. The law of displacement 
of the pole of figure assumed by Newcomb iS; however, not verified. We find 

PR : PPli ^e + e :e, 



or 



whereas Newcomb has taken 



PP':Pll=:6^€ . (a), 

PV : P'R :== E' : € (6), 



where E' denotes the ellipticity which would be induced in a sphere by centrifugal 

^ TissERANi), ' Mecanique Celeste,' vol. 2, p. 494. 

t ' Monthl J Notices,' Marcli, 1892. 

J For an account of previous Y>'ork, vide Tisserand, ' Mecanique Celeste/ vol. 2, cliap. xxx. 
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force^ if distortion were resisted by gravitation as well as by elasticity. The quantity 
e' is the eHipticity which would be induced if elasticity were the only resisting factor. 

The reason Professor Newcomb's hypothesis is at fault is that the ellipticity e, 
called by him the ^^ natural ^^ ellipticity of the spheroid, is itself partly maintained by 
centrifugal force, and that if the rotation were annulled, the spheroid would be 
distorted so that its ellipticity would no longer be e, but E, say. 

The ellipticity induced by rotation about a displaced axis has to be superposed on 
the ellipticity E about the original axis of rotation, and not on the ellipticity e. 
Thus, in (6), e should be replaced by E. Now it is obvious that e = E + E', and by 
a formula given by Thomson and Tait ^^ 



Thus 

E' 
and therefore 



1 

E' " 
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• - e 

e' 
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6 + e' 


+ 


F/ : 


Er:. 


/ 

■ € 


& %^tt 





- 9 
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With the above correction, then, the laws (a) and. (&) become identical. 

The spheroid will be distorted not only by centrifugal force about a displaced 
axis but by the relaxation of centrifugal force about the original axis. The second 
disturbing factor has been neglected by Professor New comb. 

The disturbing potential will be the difference of the rotation-potentials due to 
rotation with angular velocity co about Qz and about Oz^, that is to 

= W {^^ + ?y') - W {{^ - 'Af + (2/ + ^^i)'} ■ 

=== 0)^ {dcfcz -™ Oiyz). 

It is obvious that the equations (19), (20) are the equations for the distortion of 
an elastic sphere when distorted by forces throughout its mass derivable from this 
potential function. 

Finally the angular velocity of R about P is X, and therefore the angular velocity 
of R as viewed from P^ is 

KP' 6 ' 

w^hich is the third hypothesis made by Newcomb. 



^ ' I>ratiiral Phil osopliy,' Part II.j § 840. 
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§ 9. Conclusions, 

The existence of the Euleriaii nutation, and the fact that it would give rise to a 
variation in latitude, was first predicted tlieoretically on the assumption that the 
Earth could be regarded as a rigid body. Our present work, however, shows that 
the hypothesis of perfect rigidity, though afibrding a sufficiently close approxiniation 
to the circumstances presented by nature to specify the character of the oscillation^ 
is totally inadequate to lead to a correct determination of the period unless the 
Earth possesses a very much higher degree of rigidity than is met with in substances 
which have been subjected to exjDeriment. The only knowledge we have of the 
amount of the Earth's rigidity arises from the very vague indications furnished by 
Tidal Theory, and we must therefore have recourse to observation to determine the 
period with accuracy. This has been effected by Dr. Chandler, who, as the result 
of a discussion of a very large number of observations, has assigned 427 days as the 
true period. This period is, as the present theory requires, considerably in excess of 
the Eulerian period of 305 days. 

In § 7 we have endeavoured to obtain a numerical estimate of the effective rigidity 
of the Earth which would be consistent with this observed period, and we have 
found it to be sHghtly greater than that of steel, a result which agrees sufficiently 
closely with the requirements of Tidal Theory. Various causes however combine to 
render this result liable to a considerable amount of uncertainty. In the first place, 
our present analysis applies only to a homogeneous spheroid composed of isotropic 
material, neither of which conditions are fully realized in the case of the Earth. In 
the second place, there are probably other causes in addition to the elastic deforma- 
tions of the solid parts of the Earth, which tend to DQodify the period. In particular 
we have completely neglected the effects of the mobility of the ocean. According to 
Newcomb these effects will be small compared with the effects of elastic deformation, 
but different writers have expressed widely different opinions on the subject. Thus 
Woodward'^^ is of opinion that they alone would be sufficient to fully account for the 
observed extension of the period. Nevtcomb's view is to some extent confirmed by 
the smallness of the tide having a 427-day period which has been made the subject of 
observation by BAKHUYSENt and A. S. Christie, ;[; and as it has appeared to me to 
be quite open to question in what manner the results would be affected, I have 
thought it best not to apply any correction on this account ; the results, then, must 
only be regarded as provisional^ pending more complete mathematical investigations 
on the subject. 

In a previous paper § I have investigated the effects of an internal fluid nucleus, 

* ' Astron. Jour.,' 'Eo, 345, voL 15. 
t ' Astron. Nacli.,' No. 8261. 
% 'Astron. Jonr.,' No. 351. 
§ ' Phil. Trans.,* A, J 895. 
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and found that if the central sohd portions of a rigid spheroid were replaced by 
liquid, the theoretical estimate of the period of oscillation would be diminished. It 
might appear then that our estimate of tlie Earth's rigidity would be diminished if 
we suppose there to be a central fluid nucleus. That this is not so, I think the 
following considerations will show. 

In accordance with the results of the above-mentioned paper the effect of internal 
fluidity would be to increase the effective value of the quantity we have denoted by e. 
When the external crust is of considerable thickness, the increase in this quantity is, 
however, very slight ; thus for a crust of about 2,000 miles in thickness we find e is 
increased in the ratio 305 : 300. Now it seems that if the central solid portions of 
the Earth were replaced by fluid the increase in the value of e\ which denotes the 
ellipticity due to rotation^ would be much more rapid, and that consequently X/co 
would rapidly diminish. We conclude then that the increased effects of the elastic 
deformations would more than counteract the influence of the reduced effective 
inertia due to internal fluidity, and that with a given degree of rigidity the period of 
oscillation would be still further prolonged. The degree of rigidity of the crust 
necessary to account for a given period would also have to be increased, and as the 
estimate we have alreadv found is hip-h, the evidence furnished by the latitude- 
variation still seems opposed to the existence of an internal fluid nucleus. 

Finally we may consider the eff'ects of the Earth's viscosity. Unless this be so 
great that the present work is inapplicable^ a circumstance which seems to be quite 
precluded from the close agreement of our results with observation, the chief effect of 
internal friction w^ll be to cause the oscillation in question to gradually die out 
without producing any material change in the period. The dissipative forces arise 
entirely from the distortion of the parts of the system, and consequently no such 
forces will occur if the system be absolutely rigid throughout. Now we have seen 
that the motion consists very approximately of a rotation as a rigid body, and that 
the elastic distortions are exceedingly minute. Hence we conclude that a very small 
amount of dissipative force will be called into play, and thus if the motion is once set 
up, there appears to be no difficulty in accounting for its continuance for a very con- 
siderable period, possibly extending over several centuries. 



